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Abstract 

We analyze various models for the top quark forward-backward asymmetry (A FB ) at the Teva- 
tron, using the latest CDF measurements on different A FB s and the total cross section. The 
axigluon model in Ref. 5[ has difficulties in explaining the large rapidity dependent asymmetry 
and mass dependent asymmetry simultaneously and the parameter space relevant to A FB is ruled 
out by the latest dijet search at ATLAS. In contrast to Ref. §(], we demonstrate that the large 
parameter space in this model with a U(l)d flavor symemtry is not ruled out by flavor physics. 
The i-channel flavor- violating Z', W' and diquark models all have parameter regions that satisfy 
different Aps measurements within 1 o. However, the heavy Z 1 model which can be marginally 
consistent with the total cross section is severely constrained by the Tevatron direct search of same- 
sign top quark pair. The diquark model suffers from too large total cross section and is difficult 
to fit the tt invariant mass distribution. The electroweak precision constraints on the W model 
based on Z'-Z mixings is estimated and the result is rather weak (m^ > 450 GeV). Therefore, the 
heavy W' model seems to give the best fit for all the measurements. The W model predicts the 
ti+j signal from tW' production and is 10%-50% of SM tt at the 7 TeV LHC. Such t + j resonance 
can serve as the direct test of the W' model. 
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The prompt decay of top quark before hadronization provides opportunity to explore 
its various properties like charge, mass and spin. Given its large mass, the scale of top 
quark pair production is greater than 2m f where the perturbative QCD plays important 
role. Therefore, the top quark pair production at hadron colliders can serve as handle of 
precision test of the standard model (SM) gauge interaction, both weak interaction in its 
decay and the perturbative QCD theory of strong interaction in its production. 

From the structure of the SM, top quark is special. As a colored particle, it is the 
heaviest known particle which is copiously produced at the hadron collider. Since the top 
quark acquire its large mass through the electroweak symmetry breaking (EWSB), any of its 
properties deviated from SM would be an important signals for new physics and potentially 
indicate the origin of EWSB, which makes searching new physics in the top quark sector 
extremely interesting at both Tevatron and LHC. 

One important measurement for top quark in top quark pair production is the top forward 
backward asymmetry, which is equivalent to charge asymmetry under CP transformation 
l|. For the SM production, it involves the high precision calculation of QCD. At 0(a^), the 
bremsstrahlung amplitudes qq — > QQg carry an odd power of color charge hence have an 
odd charge conjugation parity in the interference terms among initial states radiation and 
final state radiation diagrams. There is also interference between the box diagram of 0(ot^) 
with the LO diagram that contributes to the charge asymmetry. 

CDF collaboration has recently updated the measurements on the total forward-backward 
asymmetry in top quark pair production with the semi-leptonic ti data with integrated 
luminosity of 5.3 fb _1 0. The observed total asymmetry measured in the lab frame and the 
ti rest frame are 

A FB = 0.150 ± 0.050(stat) ± 0.024(syst) (pp rest frame) 

A FB = 0.158 ± 0.072(stat) ± 0.017(syst) (ti rest frame) (1) 

which corresponds to the SM prediction based on the NLO simulation, Monte Carlo for 
FeMtobarn processes (MCFM), 0.038 ± 0.006 (in lab) and 0.058 ± 0.009 in ti rest frame 
respectively jsj]. These measurements have improved the previous results based on 3.2 fb" 1 
of A%(cos6) = 0.19 ± 0.069 and Af B (Arf) = 0.24 ± 0.014. (3] I 

1 The top quark forward-backward asymmetry has also been measured in the di-lepton channel as Afb = 

0.42 ± Q.lh(stat) ± 0.05{syst) in the ti rest frame with 5.1 fb _1 data |2|. 

2 Note that the recent DO measurement A FB = (8 ± 4(stat) ± l(syst))% is based on top-pair events that 



More importantly, with the enlarged data sample, CDF collaboration has also released 
two distributional measurements. The most interesting result is the mass dependent forward 
backward asymmetry. The mass dependent forward backward asymmetry in the tt rest frame 

A%{Mu > 450 GeV) = 0.475 ±0.112 (2) 

in comparison to the QCD correction prediction 0.088 ± 0.013. This 3.5 a deviation may 
be a strong indication for physics beyond the SM. The second measurement is the rapidity 
dependent asymmetry, which is frame independent, as 

A FB (\ Ay |> 1-0) = 0.611 ±0.256 (3) 
A FB (\ Ay |< 1.0) = 0.026 ± 0.104 ± 0.056 (4) 

in comparison to the MCFM SM prediction as A FB {\ Ay |> 1.0) = 0.123 ± 0.018 and 
A FB (\ Ay |< 1.0) = 0.039 ± 0.006. 

The ratio of the parton level asymmetries in the two different frames, which differ by 
longitudinal boost, is 

—tt = 0.95 ± 0.41 (5) 

A lt y ' 

with the error corrected for the expected correlation across frames in the NLO QCD assump- 
tion. Even though the uncertainty is still large, this close to 1 central value implies that the 
top events which contribute to the asymmetry mostly lie in the forward-backward direction 
so the asymmetries are less dependent of the longitudinal boosts along the beam direction. 
This feature is also shown in the At] dependent asymmetry Ay |> 1.0) = 0.611 which 

shows that the asymmetric events are mostly due to events with larger rapidity difference 
I Vt - m I ■ 

On the other hand, the measurement of ti cross section a t t, updated by the 4.6 fb _1 CDF 
result (with m t = 172.5 GeV), is a% p = 7.50 ± 0.31(stat) ± 0.34(syst) ± 0.15(Z theory) pb 
which is in very good agreement with SM theory prediction of = 7. 5^' j pb at NNLO []. 
Therefore, in order for new physics to generate large asymmetry without changing the total 
production cross section, the new physics contribution must interfere with the leading SM 

satisfy the experimental acceptance, which is uncorrected for effects from reconstruction or selection and 

can not be used to compare with the CDF results 4|. 



3 The latest NNLL calculation shows the a t i{m t = 173.1 GeV) = 6.30±0.19i°;23 pb fl which is significantly 
lower than the experimental results. However, we still use the old SM predictions since we do not know 
Ctf(mt = 172.5 GeV) for the latest results. 
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production of uu, dd A tt as color octet exchange in s-channel. For instance, in order for the 
s-channel massive Z' to explain the asymmetry, there is no interference between s-channel 

— Z' — — g — 

color singlet exchange uu, dd — > tt and QCD uu, dd — > tt. The asymmetry events due to 
Z' will significantly enhance the total cross section a t i at the same time and this causes a 
strong tension between fitting of Apr and a t i. This requirement implies that there are only 
two categories of candidate models to solve this anomaly. 

• First category of models contain s-channel color octet vector boson but with parity 



violation at both q — q — G and t — t — G vertices {5- 

Second category correspond to the t-channel exchange of light gauge boson of maximal 
flavor violation that couples initial state u, d quark to the third generation t quark. 



The large asymmetry can be generated via Rutherford singularity behavior 



4l6|. 



Both categories of models have their realizations in the beyond SM models. Given the 
updated measurements, especially the new distributional measurements, we discuss the cur- 
rent status of various models. In addition, the models may have other implications that 
have been or will be constrained by some direct or indirect experiments. One realization 
of the first category models is the non-universal axigluon model proposed in (5)] and it may 
receive constrain from low energy neutral meson mixings js]. However, we show that the 
flavor bound can be easily evaded by putting a horizontal flavor symmetry U{l)d- In the 
Wr models lo[ Q, since the Wr is charged under SM U(l) em , the neutral component Wr 



would inevitab 
ATLAS Dijets 



y mix with W[ and some extra U(l)x which induce a Z-Z' mixing. The new 

would severely constrain 



17] search and the Tevatron same sign dileptons 18] 
the s-channel axigluon models and the t-channel heavy Z' models. We also study the direct 
prediction at the Large Hadron Collider(LHC) using the 1 a fitting of all three asymmetry 
measurements A^ B (M t t > 450 GeV), A^ B {\ Ay\ > 1.0), A^ B (total) with the right total cross 
section. 

The paper is organized as follows. In Section HJ we presented the 1 a fitting of all 
three asymmetry measurements for s-channel color octet model (Section II A[) . t-channel 



Z' model (Section II B lj) . W' model (Section II B 21) . diquark model (Section II B 31) and the 



corresponding consequences. In Section [TT1 we calculate the production rates for the new 
particles in various different models at the Tevatron which give the bounds for those models 
and the LHC signals. In Section IHIl we consider some indirect bounds for the axigluon 



from flavor physics (Section IIII A[) and W model from electroweak precision test (EWPT) 
(Section IIII Bj) . Section ITVl contains our conclusions. 



I. UPDATED STATUS OF THE MODELS 

In this section, we discuss the updated status of the models based on the latest mea- 
surements, especially the new distributional measurements. The simulation in the following 
discussion is at parton level and leading order. The asymmetry observables are defined at 
parton level without taking into account possible reconstruction efficiency. The SM contribu- 
tion to the asymmetries from MCFM simulation have been subtracted to the corresponding 
measured values. The total cross section is obtained by multiplying a QCD fc-factor. Since 
the latest experimental value is based on m t = 172.5 GeV, for better comparison, we em- 
ploy the theory calculation at NNLO for m t = 172.5 GeV and the fc-factor is 1.3. Last, the 
differential cross section of ti invariant mass is not included as requirement in the scan since 



QCD correction 19] and cut efficiency 20] may significantly modify the shape of differential 
distribution da/d,M tr 

In the following discussion, we are mostly interested in the region where the three asym- 
metry measurements can be explained within 1 a. 

A. s-channel color octet 

The interference term between the color octet V — A gauge boson contribution and 
gluon contribution in qq annihilation pickup a term as 

J J^=^ [+ 2 9WM ] (6, 

where g s is the strong coupling, g q A is the axial component of the coupling between and 
light quarks q and g\ is that of the top quark. If the interference contribution is positive 
asymmetry, it requires that the axial coupling g\g\ < is inevitable Q. 

The s-channel models can be realized in various context. The first realization is the 
axigluon models where SU(3) C color gauge symmetry is only a remnant of SU(3)l x SU(3)r 
4 This non-universal gauge interaction potentially cause the violation of GIM mechanism thus may be 

constrained from flavor changing neutral current (FCNC) processes such as neutral meson mixings and 

we discuss its implications in the next section. 
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5|. Another realization is the models 



broken by a bi-triplet scalar and another color octet with axial coupling become massive. 
However, to achieve the g\g\ < requirement, the axigluon model has to be non- universal 
and one example is the 4-generation model proposed in 
of extra dimension theory where massive color octet Kaluza-Klein (KK) gluon couple to the 
SM quarks in the chiral form as a result of fermion profiles [21 1 . The large m t naturally 
implies that the top quark and light quarks couple to KK gluon in different way. 

One interesting feature that was discussed in the axigluon model 5| is the mass dependent 
asymmetry. Due to the opposite contribution to asymmetry between the interference term 
and new physics squared term, the asymmetry is positive when the centre-of-mass energy 
is at intermediate energy but when it is close to the threshold, the asymmetry may become 
negative. This bending-over in correlation between asymmetry A 1 and the centre-of-mass 
energy M# had been shown in the latest CDF measurements, in both the measurement with 
finite bin sizes of M t i and the measurement with below/above M t i edge. 

We use the axigluon model as one example to illustrate the feature of s-channel models 
in comparison with the updated measurement. Figure Q] shows the summary of best fit 
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FIG. 1.1a parameter region for constraints from the total asymmetry in the lab frame and the tt 
production rate; 1.5 a parameter region for the Afb in M t i > 450 and | Arj |> 1 and a(pp — > tt). 



parameter regions for total asymmetry, mass dependent asymmetry, rapidity dependent 
asymmetry, the total cross section and the last bin of da/dM t i measurements. Since the 
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total asymmetry has been reduced from the previous fitting in [5J, the 1 a region with total 
asymmetry is enlarged as shown in FigJTJ However, there is no 1 o region for the mass 
dependent asymmetry of M t j > 450 GeV or the rapidity dependent asymmetry of | At] |> 1. 
Figure [T] shows the 1.5 a parameter space for A FB {\ A77 |> 1), A FB [M t i > 450 GeV) as well 
as the total ti production rate a(pp — > ti). It is clearly shown that the axigluon model j^J 
does not consistently generate the large asymmetries in the events of M t i > 450 GeV and 
I At/ |> 1@. 



B. t-channel 

As we argued, the ratio of A pp /A tt close to one may imply that the top events are mostly 
in the forward-backward direction so the asymmetries are less dependent of the longitudinal 
boosts along the beam direction. Since the t-channel models naturally predict large number 
of events in the forward-backward region, the close to one ratio of A p v /A tt is a basic feature 
of t- channel models. 

If the asymmetry is due to new physics in t-channel physics, the interference contribution 
between new physics and SM QCD is proportional to 

C F 9 -^(u* t + sm* t + ...), (7) 
si t 

where t t = — |s(l — (3 cos 6) and 1/t expansion naturally picks up a cos 9. The t-channel 
physics naturally generates a large asymmetry in the ti system. In addition, the maximal 
asymmetry is generate at the Rutherford singularity where 8 = which corresponds to very 
high centre-of-mass energy. One would then expect the positive correlation between A FB 
and Mu. 

The t-channel Z' model in 0, [12I proposed a color singlet neutral gauge boson with max- 
imal flavor violation between first and third generations and the new contribution interferes 
with the SM uu — > ti. Similar to the Z' model, instead of neutral current exchange in 
t-channel, there is also a proposal using charge current exchange in t-channel as flavor vio- 



lation W. The interference effect is reduced since it's only the dd initial state [uj 



3. Such 



5 One can use the general color octet vector boson with V — A interaction to fit the two distributional 
asymmetries and total asymmetry [3|. The general results together with the most recent bounds in the 
dijet channel from ATLAS will be presented elsewhere. 
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lavor violation gauge interactions may be realized in horizontal gauge interaction models 



131 ] for neutral current or generalized left-right model 14| for charged current. 



A Higgs-like scalar with maximal flavor violation 22( would generate a large negative 
asymmetry due to the helicity-flip in the Yukawa coupling. The spin conservation in the 
= direction requires the top quark to move backward. To resolve this, the fermion- 
number violating diquark scalars with maximal flavor violation was proposed ll|, [l5|, |23 |. 
Diquark scalar can be 3 ® 3 = 6 © 3 under SU(3) C and has fermion-number violating 
coupling as t c u<fi or t c d<fi. Such diquark scalars with flavor violation can also be realized 
in various BSM contexts, partial unification models or supersymmetry. For instance, im- 
parity violation supersymmetric standard model which contains the baryon number violating 
coupling, e a p 1 u c a d c pdfc t [24j, the down type squark di can mediated w-channel dd — > tt that 
interferes with the QCD dd — > tt. 

All the three proposals can in principle predict large positive asymmetry in tt production. 
In the following paragraphs, we examine the numerics to see whether the models can explain 
the three asymmetry measurements and the total cross section at the same time. 



1. Z' 



We first examine the first proposed t-channel model, Z' 9]. To minimize the constraints 
from low energy, the authors proposed a right-handed coupled Z' with large coupling between 
u and t. The parameter region for 1 a fitting of all three asymmetry measurements as well 
as the total cross section for light Z' mostly below t threshold is presented in Figure El Due 
to large destructive interference, the total cross section is always smaller than the measured 
value. This result is also shown in the NLO calculation of Z' model [if]]. The best fit points 
for heavy Z' by requiring 1 a fitting for all the three asymmetry measurements are listed 
in the Table [B The corresponding tt cross section are also below the 1 o total cross section 
and the best points are towards heavy masses of (9(700 GeVV 



One more complication which has been discussed in [9-|ll| is that the events in the 
i-channel exchange tend to be in high energy region which significantly increase the tail 
of da/dMu, especially the last bin (800 GeV-1.4 TeV) in da/dM t r. The QCD correction 
may change the shape and lower the contribution in high energy 19j. In addition, the t- 
channel kinematics implies that the top quark events at high energy are mostly in the larger 
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FIG. 2. Parameter space scan of 1 a for all three asymmetries, A 1 ^ 1 , A F b(\ A77 |> 1) and 
A FB (Mu > 450 GeV) is shown in red. The 1 a fitting for A FB {M t i > 450 GeV) and the total 
cross section a t j. 
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24.7% 


6.9 


750, 1.9 


15.9% 


34.7% 


23.2% 


6.9 


CDF 


5.7%-16.7% 


27.5%-50.0% 


23.1%- 74.5% 


7.5 ± 0.48 



TABLE I. 1 a benchmark points for all three asymmetry measurements, /c-factor = 1.3, m(=172.5 
GeV for cr t j. 



rapidity region while the selection cut are more efficient for the central events. Consequently, 
the cut efficiency at high invariant mass is quite low 20|, which may further decrease the 
effective total cross section. Polarization of top quark in the events sample also effect the 
cut efficiency. 



2. W 



To resolve the tension between cross section and total asymmetry in the Z' model, the 
charged current process in t-channel may give better fit which has smaller interference effect 
due to the dd initial^state. We plot the allowed parameter regions for the t-channel charged 
current model in 
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FIG. 3. Parameter space scan of 1 a for all three asymmetries, A^g 1 , Afb(\ A?] |> 1) and 
Afb(Mu > 450 GeV) is shown in red. The 1 a fitting for A FB (M a > 450 GeV) and the total 
cross section <r t j. 



The 1 a asymmetry region of all three measurements corresponds to a larger total cross 
section which is outside the 1 a fit of latest a t t measurement. However, various efficiency 
effects discussed in the last paragraph of Z' session may significantly reduce the measured 
cross section. 



3. Diquark 

We use the anti-triplet diquark that couples to t c u<j) to illustrate the feature. Similar to 
the W case, there also exist diquark scalars whose couplings are of t c d(p and these diquark 
scalars contribute to dd —> ti instead. 
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Figure H] gives the 1 a fitting for the anti-triplet diquark scalar with maximal flavor vio- 
lation. The 1 a region also exists for the anti-triplet diquark scalar for all the measurements 

8 
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FIG. 4. Parameter space scan of 1 a for all three asymmetries, A 1 ^ 1 , Afb(\ Ar/ |> 1) and 
A FB (Mu > 450 GeV) is shown in red. The 1 a fitting for A FB {M a > 450 GeV) and the total 
cross section <r t j. 



in asymmetries A FB . But the corresponding total cross section a t i are also larger than the 
measured value over 1 a . In addition, the best-fit for cross section and the mass dependent 



asymmetr y is over 1 TeV which makes t 
shown in ll|. The latest simulation by 



re dajdMti measurement very difficult to fit as 



20[ also showed that the tt events generated by 



diquark scalar had a higher cut efficiency at high energy therefore the anti-triplet diquark 
fitting is worse than the W. 



II. IMPLICATIONS AT THE TEVATRON AND LHC 

After fitting the top forward backward asymmetries in different kinematical regions, we 
discuss the other Tevatron bounds for the models and the LHC predictions that can be soon 
tested in this section. 

The Large Hadron Collider (LHC) is a proton-proton collider with centre-of-mass energy 
7 TeV in the first two years running. Unlike at Tevatron where the axigluon effect only 
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appears as interference. The color octet axigluon of 0(1 TeV) can be directly produced 
at the LHC and decay into dijet or tt. With significant decay branching ratio (BR) to ti, 
it provides additional handle to search it. The study of axigluon at the LHC has been 
performed by 7|. ATLAS collaboration has recently released the search for dijet resonance. 
The latest data has ruled out axigluon from 0.6-2.1 TeV by assuming axigluon coupling is 
only g s . The axigluon model in |5j has a even larger coupling comparing with the ATLAS 
paper and therefore, the model receive much more server constraint. 

For neutral gauge boson like Z', the flavor violating vertex of ut will lead to large uu — > tt 
or uu — > tt scattering with Z' exchange in the t/w-channel. The same-sign positive top quark 
pair (uu — > tt) becomes particular interesting at the LHC given its large w-valence quark 
parton flux [251 ] . 

In addition, with large ut coupling, the tZ' associate production is not negligible. Since 
Z' equally decays into ut and tu, the associated production tZ' or iZ' will contribute to 
tt + j, tt + j and ti+j final states. Again, since the LHC is proton proton collider, the tt + j 
dominates the same-sign top production. The ti + j will appear in the inclusive ti search. 
Since the 1 a parameter space of all the asymmetry constraints corresponds to smaller ti 
pair production, the additional ti + j may in principle help to ease the tension at Tevatron. 
However, if it significantly contribute to the ti, the same amount of same-sign top quark 
will arise. 

Figure [5] (a) gives the pp — > tt production rate at Tevatron and the 7 TeV LHC with the 
Z' in the 1 a fitting for all three asymmetry measurements. The pp — > tt + tt at Tevatron 
is between 0.7-1 pb for these best fit points. CDF measured only 3 events for 2 fb _1 [l8| 
with the acceptance range from 1.5% to 3%. The best fit points all predict 15-30 same-sign 
pure leptonic top events before selection cut but with one b-tagging. Even though these 



events from t-channel vector boson exchange may suffer fromji 
to the t-channel light scalar exchange considered in Ref. 



ow cut efficiency comparing 
, the Z' model is strongly 



constrained by the same-sign top quark scattering. At Tevatron, the same-sign top due to 
tZ' associate production is then much suppressed at the Tevatron due to significant phase 
space suppression. 

The uu — > tt scattering get significantly enhanced at the proton-proton collider LHC. 
The production rate can reach 200 pb. Therefore, even at very early running of LHC with 
about 30 pb -1 data and requiring two b-tagging jet, the event number before kinematic cut 
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FIG. 5. (a) a(pp -)• tt) at Tevatron and the 7 TeV LHC; (b) a(pp ->• tZ' + t + Z') at Tevatron 
and the 7 TeV LHC. Both (a) and (b) are based on model parameter of Z' in the 1 a fitting of all 
three asymmetry measurements as listed in Tabled! 



is about 70 and the same-sign top quark tt events is expected to be (9(10). 

For W or diquark scalars with flavor violation, since they are electrically charged, it will 
only contribute to tt as at Tevatron. However, since the W or diquark <fi has a large dt or 
ut coupling, the dg — > tW' or ug — > i<p production is significant as shown in 10, O, l^ . 
With W and diquark scalars of typically above top quark threshold, they can decay into t 
plus one hard jet. The signal is then tt plus one hard jet and should appear in the inclusive 
tt searches. The diquark case has already been calculated in our early paper 11]. 

Figure [6] gives the production of top quark plus W at hadron colliders. For Mw' lighter 
than 400 GeV, the production rate is about 0.1-1 pb at the Tevatron. As we discussed 
in the previous section, the 1 a fitting parameter space for all the asymmetries constraint 
corresponds to the larger cross section region. |f The new contribution to tt + j will increase 
the tension between Afb and a t t- For heavy W above 400 GeV, due to large phase space 
suppression, the production rate at Tevatron can then be neglected. However, at the LHC, 
even with 7 TeV centre-of-mass energy, the production rate is O(10 pb). Even though 



261 ] claims that the single top production at Tevatron puts a strong constraint on the W' — u — b coupling. 



However this constraint does not apply to general W' models. 
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FIG. 6. pp -»• tW ->j + tt 

gg — > ti dominates the ti at the LHC which does not interfere with W, the tW itself may 
significantly increase the ti rate. 

III. INDIRECT CONSTRAINTS FOR MODELS 
A. Axigluon with flavor protection 

As shown in the previous section, only non-universal axigluon models can provide the 
positive asymmetry. Being a color octet with strong coupling strength, this GIM violation 
axigluon will then lead to significant flavor changing neutral current (FCNC) effect. 



£ = igtf-f (HI) P L qiT a G; + ig R i?<f (H%) P R q j T a Gl 



(8) 



Flavor violation thus can arise from the non-universal gauge couplings due to the rotation 
between mass eigenstate and gauge eigenstate. 



u L = V£u L , d L = V£d L , u R = V^u Rl d R = V^d 



R 



The effective coupling in horizontal space is 



V^H U L V^ V^H-VS, V^HiVh V^H d R Vi, 



(9) 



(10) 



14 



The rotation from mass eigenstate to gauge eigenstate for up and down type quark 
respectively is completely unmeasurable in the weak interaction. The only observable is the 
mixing in charge current transition which is categorized as CKM matrix 

To avoid flavor violation in the down sector, one may introduce a U(X)d symmetry 28] 
which acts only on the down sector with different eigenvalues for different generations, but 
does not distinguish the handedness of the quarks. Then the down-quark sector are diagonal 
with VI — Vft — 1 so that there is no FCNC at all in the B s , or neutral K system. For 
simplicity, we take further the rotation matrix of right-handed up-quark sector as V£ = 1. 
Then one can explicitly determine the left-handed up-type quark rotation based on the 
known CKM matrix using V£V^ = Vckm 

Vt = Vckm (11) 



Nowadays the Wolfenstein parametrization 29[ is widely used to express the CKM matrix 



in terms of four parameters (A, A, p and 77). To keep the unitarity of CKM matrix to all 
orders of A, we adopt in the following a definition of Wolfenstein parameters proposed in 
[30I . Then the effective coupling between up quark and charm quark is 

{VtfHlVZU = ~A 2 X 5 ^v ~ P + 1) • (12) 

Under the assumption of above rotations, the FCNC operators only arise in left-handed 
and mixing between first and second generation up-type quarks as 

- ^i-f^M-rA) + , (13) 

where the following decomposition satisfied by the color SU(3) fundamental representation 
has been implemented 



Under Fierz transformation 



K^e = o^eVr - JaAe • (14) 



(utr4M^C a L ) = (UI^CI)(U^4) , (15) 

the effective AC = 2 Hamiltonian can be expressed as 

^ A r=c(/i)(«M)(«M (i6) 

7 It was argued in [|| that the axigluon model in 5 suffers serve bounds from Bd mixing. However, the 
calculation seems to be done by assuming both up and down quark sectors transform like CKM rotation. 
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and the leading order Wilson coefficient at the scale ma is 

C(mo) = " (17) 

It means that the D° — D° mixing in this axigluon model has A 10 suppression due to CKM 
rotation. The RGE running of the above Wilson coefficient is well knowrjf , 

Notice that the hadronic matrix element of AC = 2 operator is 

(D°\uirc a L uhA\D°) = \f D m 2 D B D {p c ) . (19) 

Just like the B° - B° mixing case, one may define the renormalization group invariant 
parameter Bp by 

B D = (a s (fi c ))- 6 / 23 B D (fi c ) , (20) 

which should be 0(1). Then the axigluon induced AC = 2 effective operator contributes to 
mass difference of neutral D system as 



m 8irf*m D B D a a {m G )A*\ 10 ((l - pf + ^ 
9 th3 



Am D = a^r -'V — — ' ^ " ! •> > . (21) 



The axigluon model can also induce AC = 1 effective operator which would in principle 
affect D° — D° mixing by 

M--F) = —(D \H5 c f = 2 \D°) + — s P - ' e// ' A ' e// ' -. (22) 

2 i„ 2m D y 1 e// 1 1 2m D ^ m D -E n + ie K J 



Actually the experimental observation of comparably large mass and width differences 31] 

x= ^ = 0.98ig£% , y = ^ = (0.83 ± 0.16)% (23) 

strongly implies that they are dominated by the long distance effects of the SM AC = 1 
operators. Therefore the axigluon induced AC = 1 terms could be safely neglected as they 
should be much smaller than the tree-level SM AC = 1 terms. 
Taking the Wolfenstein parameters as 32] 



A = 0.812 , A = 0.2254 , p = 0.148 , r] = 0.351 (24) 



The RGE running is actually dependent on Nf, the number of active flavor via = 11 — 27V//3. From 
the scale tug down to /i c , Nf changes correspondingly from 6 to 4. But numerically this effect is small 
and we will simply take Nf = 5 in the RGE running. 
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and f D = 207 MeV [3l|, we obtain 

(^£) = 0.082% (i-^YY^) (25) 

which is roughly one order of magnitude smaller than the experimental result. 

B. Electroweak constrains on the W' model 

In general, the W must generate its mass through gauge symmetry breaking, then some 
other neutral component in the W symmetry breaking sector (for instance W R in the SU(2) R 
symmetry breaking) would inevitably mix with and some extra U(l) so that the W r± 
would be charged under U(l) em . As a consequence, there is a large Z-Z' mixing which is 
constrained by the electroweak precision test. In general, the bound from EWPT is subtle 
since different fermions W/Z boson couplings are modified in different ways which may even 
depends on models and a careful global fit is needed. The full results for the W model to 
explan the top forward backward asymmetry will be presented elsewhere. Since the overall 
modification for fermion Z boson coupling is small (except for some right-handed quarks 
charged under SU(2) R ), we only consider the tree-level Z-Z' mixing as a rough estimation. 
For observables that strongly depends on u/d/bR-Z coupling, such as g R , Qw(Cs), etc., their 
deviations from the SM results are still at the same level as Z mass which is transmitted 
into W boson mass. 

We can start to consider a simple SU(2) R x U(l)x x SU{2) L model to estimate how 
large is the electroweak constraint for the Z' — Z mixing. The SU(2) R is separated from 
SU(2)l to avoid the troublesome W — W mixing. The two double Higgs field ht and 
are charged under SU(2)l x U(1)x and U(l)x x SU(2)r respectively. The higgs fields 
get their vacuum expectation values (hi) = ul and (h R ) = ur which spontaneously break 
SU(2) R x U(l)x x SU(2)i into the diagonal group U(l) em . In order to raise the Z' mass 
so we have less constrain from the Z-Z' mixing, we choose hn transform as a triplet under 
SU(2)r so mz> = \f2mw 1 - The gauge quantum number for h L and h R are (0, 1, 1/2) and 
(1,2,0) under SU(2) R x U(l)x x SU(2) L respectively. For SM fermions, at least the quark 
doublet (t,d)R is charged under SU(2) R (It is possible to have some extra hidden fermions 
charged under SU(2)r and U(l)x to cancel the gauge anomaly). For the rest SM fermions, 
their quantum number is the same as the SM one if one replace their hyper charge as the 
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U(l)x charge. The quantum number for (t, d) R and (u, b) R under SU (2) R x U{1) x x (2) l 
are (1/2, 1/3, 0). 

The kinetic term for the link fields Tr[(D At h L )^(D At h L )] + Tr[(D At h R )'''(D /J h R )] becomes the 
mass terms for the massive gauge bosons. The mass matrix of the gauge bosons is 

1 2gW p -2a»ovui o \ {a*\ 



-2g R g x u R 







-2g R g x u\ g 2 x (2u 2 R + u 2 L ) -g X 9LU 2 L 



\ 



(26) 



-g x g L u z L 

We introduce the parameter e = u 2 L /2u R <C 1 which shows that the right-handed symmetry 



breaking is only a perturbation. 

This matrix can be diagonalized by means of an orthogonal matrix which we shall call 

R: 







(a \ 










(27) 


\ A i) 









where the mass eigenstates are denoted by A, Z, and Z' . The eigenstate A is massless and 
identified as the photon. The couplings of our theory are related to the electric charge by 

e e e 



9R = ^— , 9x = ^— , 9l = . a (28) 

sin <p cos u w cos (p cos Uw sm t) w 

where 9w is the weak mixing angle (in the limit e — > 0) and is an additional mixing angle. 

The other two eigenmasses are 



rrir 



rrir 



u 2 l (9y + 9 2 l) 



1 
2 

2/2 i 2i 

2 n i?(^ + &o 



+ 9 2 x ) 2 



-u 2 l(9y + 9 2 l) [l-"in 4 



(29) 



(9 2 R + 9 2 X ) 2 



2«r(9r + 9x) [l + esin 4 0] , (30) 



where we have dropped 0(e 2 ) terms and 



9 2 y 



9l 



9\ 



(31) 



Clearly, Z is identified with the SM Z boson while Z' is referred to as the heavy Z boson. 
For small e, the mixing matrix R has the following approximate form: 

/ c^A^^a..- „™ a... ^ a... \ 



R 



sin <p cos Ow 

„3 . 



i • /i , sin cos 

sin <p sin t>w + e 
y — cos <fi + 



sin % 

e cos sin 4 r£ 



cos <p cos f 
cos </> sin 6*w/ — e 



cos sin 



e cos 2 d> sin 3 i 



Sill C/v7 

— cos 8w 
-e cot 0jy cos (f> sin 13 0y 



(32) 
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FIG. 7. The excluded region on the parameter spaces of Wr model based on tree level contribution 
to T parameter alone from the new physics at 90% C. L. The green plus yellow region are is the 
excluded region for SM Higgs mass = 114GeV while the red plus yellow region is the excluded 
region for SM Higgs mass = 150GeV. 



from which it is simple to derive the SM fermion couplings. The SM fermion and Higgs 
couplings to Z and Z' can be written as 

gR sin sin W T R + g x cos sin OwQx — 9l cos O w Tl 

(sin 2 9 W Q -T'l + e sin 2 cos 2 <f)T R - e sin 4 <pQ x ) , (33) 



sin 9w cos 9\y 



gR(— cos0 + ecos0sin (j))T R + gx(sin0 + ecos 0sin 4>)Qx — 9l(— ecot6w cos0sin 0)T L 



cos ^ 3 sin e cos sin 3 ^ 3 



-T R + —T-Q x + ,\ y (-T 3 + Q sin 2 M . (34) 
cosfcV \ sin<p cos0 sin uw J 

In the limit of large SU(2)r breaking vev (e 1) and small mixings (0 — > 0), the Higgs 
current can be approximated as (we drop the e cos 2 sin 3 term) 

J%,(h) = --g x sin ^D^h) + /i.e. , (35) 
which induce a dimension six operator 

a h O h = sin4 <K^^) 2 • (36) 
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which coincident with Eq. ( 129]) that Am| = — esin0 4 m|. We can calculate the correspond- 
ing T parameter from the tree level gauge boson mixing, 

rp = _a>y = e sin 4 <p 
a f a f 

Using the SM model mz, Gf (the life time of r) and a = e 2 /4ir as the basic input param- 
eter, we can calculate the allowed parameter space including the Higgs radiative corrections 
according to the most recent results: S = 0.03 ± 0.09 and T = 0.07 ± 0.08 (with 87% strong 



correlation) 



3l| . The results are presented in Fig [7J We can see that for sufficient heavy Z' 



and strong coupling gn (for instance, gn = 2, mz> = 900 GeV which is used in Ref. [20j), it 
is well above the excluded region. 

IV. CONCLUSIONS 

We discuss the feature of various models for the top quark forward-backward asymmetry 
anomaly at Tevatron, using the latest CDF measurements on total asymmetry in lab frame 
ApB lab , the rapidity dependent asymmetry Aps{\ A?7 |> 1), the mass depedent asymmetry 



A FB (Ma > 450 GeV) and 
The axigluon model in 

and the mass dependent asymmetry simultaneously. In addition, the latest dijet search 17] 



the total tt production cross section o t i- 

5]] has difficulty in explain the large rapidity dependent asymmetry 



at ATLAS has ruled out the parameter region that is relevant to top Afb- On the other 



hand, in contrast to the conclusion in Ref. the model itself does not suffer from the 
flavor Bd mixing under flavor protection U(l)d and a careful calculation shows that their to 
D° — D° is still one order lower than the current experimental bound. 



The t-channel Z' 9|, W 1Q[ and anti-triplet diquark 11] models all have parameter re- 
gions that satisfy all three asymmetry measurements within 1 a. However, the corresponding 
production cross section predicted by the 1 a asymmetry requirement in the Z 1 model are 
always significantly below the 1 a of cross section measurement. The best fit point of Z' 
is about 700 GeV with purely righthanded coupling g^ t ~ 1.8 which corresponds to 6.9 pb. 
However, this best fit point will generate a large number of same-sign top quark events at 
Tevatron which is at least five times larger than the SM prediction. We conclude that the 
Z' model is very difficult to be consistent with all the measurements. 

Both W and anti-triplet diquark models predict the cross sections are larger than the 
measurement but various factors can lower the survival efficiency after cuts in these models 
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to ease the tension between asymmetry and cross section. The best fit point for anti-triplet 



diquark lies in very high mass region and with better survival efficiency [20] , it is difficult to 
fit the differential cross section da/dM^. A rough estimation for W model shows that the 
bounds from electroweak precision tests are weak due the heavy Z' and strongly coupling 
g^. Therefore, we conclude that the best model is the t-channel W model at the current 
stage. To test such model directly, we also use the 1 a asymmetry parameters to compute 
the production rate of ti + j from tW at 7 TeV LHC and the production rate is 10%-50% 
of SM ti. 

Last, we want to mention that the latest NNLL calculation a t i{m t = 173.1 GeV) = 
6.30 ± 0.19^Q23 pb []J is significantly lower than the experimental results. If the result does 
not significantly change for m t = 172.5 GeV which is used for Tevatron experiments, then 
the fits for t-channel W and anti-triplet would be better while the t-channel Z' would be 
worse. 
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NOTES ADDED 



While this work was being delayed by the huge earthquake in Japan, Ref. 20[ appeared, 
which overlap with ours in the study of fitting different models on different top forward 
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backward asymmetries. Our results agree with quantitatively with theirs for fitting different 
top forward backward asymmetries and the total ti cross section. However, we notice that 
the axigluon model and the heavy Z' model are severely constrained by the dijet search at 
the LHC (ATLAS) and the same sign dilepton search at the Tevatron (CDF). Therefore, we 
conclude that the heavy W model is the most promising one at present. We also consider 
the indirect bounds for different models. 
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